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I review the status of fixed order jet Monte Carlos and briefiy discuss the prospects for 
next-to-next-to-leading order calculations. I present a general purpose next-to-leading order 
Monte Carlo program for four jet event shape observables in electron-positron annihilation and 
present some estimates of the light jet mass and narrow jet broadening distributions. Finally, 
I discuss an estimate of the strong coupling constant using the measured energy evolution of 
the average value of event shape variables such as thrust and heavy jet mass. 

1 Introduction and status of fixed order parton level calculations 

In the last decade there has been enormous progress in using perturbative QCD to predict and 
describe events containing jets. At the simplest lowest-order (LO) level, each jet is the footprint 
of a hard, well separated parton produced in the event. Although the predicted rate is sensitive 
to the choices of renormalization and factorization scales, qualitative comparisons of data and 
theory are generally very good. Techniques for computing tree level multiparton scattering 
amplitudes are very advanced and reliable numerical Monte Carlo programs exist (see Table |l]) 
for processes such a.s pp W + 4 jets which provides one of the main backgrounds for top 
quark detection at the TEVATRON. However, the bulk of events populate topologies with fewer 
jets and here a more quantitative description is required. To date, this has been largely achieved 
by improving the theoretical prediction to next-to- leading order (NLO). 

The addition of NLO effects produces four important improvements over a LO estimate. 
First, the dependence on the unphysical renormalization and factrorization scales is reduced 
so that the normalization is more certain. Second, we begin to reconstruct the parton shower 
so that two partons may combine to form a single jet. As a result, jet cross sections become 
sensitive to the details of the jet finding algorithm, particularly the way in which the hadrons 
are combined to form the jet axis and energy, and to the size of a jet cone. This sensitivity 
is also seen in experimental results. Third, the calculation becomes more sensitive to detector 



limitations because radiation outside the detector is simulated. This can change leading order 
results considerably for quantities such as the missing transverse energy in events containing 
a W boson. Fourth, the presence of infrared logarithms is clearly seen and regions where 
resummations are needed to improve the perturbative prediction can be identified. 

Where possible, NLO estimates of QCD cross sections have become de rigeur and have 
allowed many tests of the underlying dynamics and measurements of, for example, the running 
strong coupling constant. For all types of experiment, NLO Monte Carlos are playing a vital 
role in making detailed comparisons with hadronic events. Some examples of NLO programs 
are listed in Table |l[ Recent examples where disagreement between theory and experiment have 
yielded great excitement are the single jet inclusive data at CDF and the + 1 jet/ W +.0 jet 
ratio measured by DO. For updates on these issues, I refer you to the talks by Elvira □ and 
Summers 0. 

However, one loop matrix elements rapidly become more difficult to compute as the number 
of external particles increases. In the last few years, a number of problems connected with 



and the one-loop scattering 
a have been computed. These 



one-loop integrals with five external particles have been solver 
amplitudes for five partonQ and four partons and a vector boson E 
matrix elements are relevant for pp ^ 3 jets, e~^e~ — > 4 jets, ep ^ 3 + 1 jets and pp ^ V + 2 jets. 
This latter process is likely to be very relevant at the TEVATRON in Run II, where associated 
Higgs production pp ^ V + H followed hy H ^ bh is an eagerly anticipated discovery channel 
for the Higgs boson. NLO Monte Carlo programs for some of these processes are starting to 
appear Bi andi0 (see Table |^) and comparisons with experimental data are underway. 

However, while the NLO program has been extremely successful in reducing the theoretical 
uncertainty, the improvements in the data are even more impressive, to the extent that the 
theoretical error tends to dominate. For example, in extracting the strong coupling from event 
shapes at LEP, the renormalization scale uncertainty engenders approximately a 5% uncertainty 
in as{Mz) while the experimental error is typically 2%. One way to improve the theoretical 
predictions is to incorporate next-to-next-to-leading order (NNLO) effects. This would reduce 
the renormalization scale uncertainty and also allow a better interface between the theoretical 
and experimental jet finding algorithms - now three partons will be able to merge to form the 
jet. Although a complete NNLO jet calculation is some way off, and for 2^2 processes the 
two-loop double box integrals are not even known yet, some encouraging steps have been taken 
in this direction. For example, the singularity structures when two partons are unresolved has 
been examinedS. More impressive is the recent calculation of the two-loop four gluon scattering 
amplitudes in a = 4 supersymmetric version of QCD0. While this result may not be directly 
relevant, the techniques developed to obtain it undoubtedly will be and I expect rapid progress 
in the evaluation of two loop QCD matrix elements. Much work will still be needed to turn the 
amplitudes into cross sections, but I predict that numerical NNLO programs will be available 
for pp ^< 2 jets and e~^e~ — > 3 jets will be available within the next five years. 

In the remainder of this talk, I will focus on two recent examples of theoretical progress in 
jet physics. First, I will discuss some new results for four jet event shape observables using the 
recent one loop matrix element calculations for 1/ — > 4 partons i'i. Second, I will discuss how 
the theoretical error on the strong coupling may be reduced by considering average values of 
event shapes at a range of different centre-of-mass energies. 



2 4 jet observables 



As mentioned earlier, the one-loop matrix elements appropriate for the decay of a vector boson 
into four partons have recently been computed by two groups BQ. These matrix elements have 
been combined with the tree level amplitudes for the decay of a vector boson into five partons 
to compute the NLO corrections to a variety of four jet rates and event shapes by two groups; 
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Table 1: Available multi parton scattering matrix elements and some commonly used fixed order Monte Carlo 

implementations. 



Algorithm 


ycut 


MENLO PARC 


DEBRECEN 


EERAD2 


Durham 


0.005 
0.01 
0.03 


(1.04 ±0.02) • 10-1 
(4.70 ± 0.06) • 10-2 
(6.82 ±0.08) • 10-3 


(1.05 ±0.01) • 10-1 
(4.66 ± 0.02) • 10-2 
(6.87 ±0.04) • 10-3 


(1.05 ±0.01) • 10-1 
(4.65 ±0.02) • 10-2 
(6.86 ± 0.03) • 10-3 


Geneva 


0.02 
0.03 
0.05 


(2.56 ±0.06) • 10-1 
(1.71 ±0.03) • 10-1 
(8.58 ±0.15) • 10-2 


(2.63 ±0.06) • 10-1 
(1.75 ±0.03) • 10-1 
(8.37 ±0.12) • 10-2 


(2.61 ±0.05) • 10-1 
(1.72 ±0.03) • 10-1 
(8.50 ±0.06) • 10-2 


JADE-EO 


0.005 
0.01 
0.03 


(3.79 ±0.08) • 10-1 
(1.88 ±0.03) • 10-1 
(3.46 ± 0.05) • 10-2 


(3.88 ± 0.07) • 10-1 
(1.92 ±0.01) • 10-1 
(3.37 ±0.01) • 10-2 


(3.87 ±0.03) • 10-1 
(1.93 ±0.01) • 10-1 
(3.35 ±0.01) • 10-2 



Table 2: The four-jet fraction as calculated by MENLO PARC, DEBRECEN and EERAD2, for the different 
jet recombination schemes and varying i/cut- The rate is normalized by the 0{a.s) total hadronic cross-section, 

CTtot = CTQ (1 + Qs/tt). 



the programs MENLO PARC by Dixon and Signer i and DEBRECEN by Trocsanyi and Nagy 
liy. Here we report on results obtained using a third numerical implementation of these matrix 
elements to compute infrared safe four jet observables, EERAD2 H. This program uses the 
'squared' one-loop matrix elements of i together with squared tree level matrix elements for 
7* — > 5 partons. Both four and five parton contributions are singular in the infrared limit and 
there are several ways ofnerforming the cancellation B3'C3'EjIH3. We use a modified version of the 
generic slicing approachtJ that does not depend on the slicing parameter ymin- 

As a check of the numerical results, Table ^ shows the predictions for each of the three 
Monte Carlo programs for the four jet rate for three jet clustering algorithms; the Jade-EO, 
Durhamc^, and GenevaE3 algorithms. We show results with agiMz) = 0.118 for three values of 
the jet resolution parameter ycut- There is good agreement with the results from the other two 
calculations. 

In addition to the jet rates, there are also NLO predictions in the literature for four-jet 
event shapes such as the D parameter, Acoplanarity and the Fox- Wolfram moments 0. These 
variables vanish as the three jet limit approaches; i.e. as the event becomes more planar. Of 
course there are many other event shape variables that vanish as the event becomes more three 
jet like. Here, I focus on observables derived by dividing the event into two hemispheres H\ and 




Figure 1: The (a) light hemisphere mass distribution l/a^adda /d{M'f^/ s) and (b) narrow jet broadening distribu- 
tion l/ahadda/dBjnin at LEP energies. The LO (NLO) prediction is shown dashed (solid) for fi = Mp-. The LO 
(and NLO) prediction using the FAC scale is shown as a dotted line. The data is taken fromEj. 



H2 according to the orientation of the thrust axis n defined by, 

r = max%%^. (1) 

Efe \Pk\ 

Particles that satisfy pi.n > are assigned to hemisphere Hi, while all other particles are in H2. 
The light hemisphere mass, m\/s, and the narrow jet broadening, i?min) are defined by, 

2 

mi 




B^ir, = min ^'^^ ' ,1 , , (2) 



are both non-zero when there are at least four particles in the final state. These variables are 
related to the heavy jet mass m?^/s and the wide jet broadening -Bmax which are obtained by 
maximising the quantities in Eq. |2|, and which require at least three particles in the event to be 
non-zero. 

The LO and NLO predictions for m\/s and i?min with a renormalization scale ^ = Mz 
and as{Mz) = 0.1194 are shown in fig. 0, together with data from §1. At small values of 
the observable, we see the presence of large logarithms which must be resummed to obtain a 
meaningful result. At larger values, we see that the next-to-leading order corrections are large 

tnd approximately 100%. Similar large effects have been noted for other four jet observables 
'113. While such large effects may make the perturbation theory appear unreliable, we note that 
effects of similar size have been noted for three jet event shapes like thrust This may be a 
sign that large higher order or power corrections are present. Or it may simply indicate that the 
physical scale /U = Mz is a very poor choice of the renormalization scale, which itself engenders 
large ultraviolet logarithms in the higher order terms. In fact, such ultraviolet logarithms may be 
resummed using the FAC or fastest apparent convergence scaleS. This scale is usually much less 
than the physical scale (and particularly so when the next-to-leading order corrections appear 
large) and may be disfavoured for that reason. However, the prediction using this scale (shown 
as dotted lines in fig. ||.) does lie much closer to the data. 



3 Strong coupling from the energy evolution of event shapes 



A vital ingredient in perturbative predictions is an accurate knowledge of the strong coupling 
constant. This can be determined via analysis of event shape variables at LEP. Consider the 
observable R{Q) with a perturbation series and leading power correction that describes the 
hadronization phase of the hard scattering of the form, 

R{Q) = a + ria^ + rgo^ + . . . + ^(1 + Aia + . . .), (3) 

where a = as{n)/TT denotes the renormalization group improved coupling. Note that the nor- 
malization is simply such that the perturbative expansion begins with unit coefficient. An 
example of such a variable is R{Q) =<1 — T> /1. 05, where in the MS scheme with ^ = Q 
and Nf = 5 active quark flavours the NLO coeflicient is ri = 9.7000. The NNLO coefficient 
r2 is as yet unknown. The precise form of the power corrections is as yet not fully understood, 
but, for the purposes of comparison with data, may be parameterized in a variety of different 
waysc2l. To extract as{Mz) from the data, we just truncate the perturbative series for a given 
renormalization scale fi = xQ (which is typically x = 1). In other words, we assume the higher 
order term ^2 = 0, = etc. Then, by comparing with experimental data, we solve for a. A 
recent analysis for <1 — T> using the power correction model ofS finds, 

as{Mz) = 0.1204 ± 0.0013 t'oZfo ^-omll (4) 

(with a x^/d.o.f of 42.6/24) where the first error is purely experimental. The second and third 
errors come from varying the theoretical input parameters, first varying the renormalization 
scale X between 0.5 and 2 and second the parameters of the power correction model. Clearly the 
estimate of the theoretical error is dominated by the renormalization scale uncertainty. Similar 
results are presented inB3. 

Alternatively, we may avoid the renormalization scale errtirely and directly write an expres- 
sion for the running of R{Q) with Q in terms of R{Q) itself E3'E3, 

JTj> 

-6i?2(l + cR + p2R^ + ...) + KoR~^/^e~^/^^il + K^R + ...) + ... 



d\ogQ 

= -bp{R). (5) 



Here b and c are the first two universal terms of the QCD beta-function, 

33 - 2iV/ _ 153 - 19Nf 

6 ' 126 ■ 

The quantity. 



(6) 



P2=r2 + C2-ric-rj, (7) 

is a renormalization scheme and renormalization scale (RS) invariant combination of the NLO 
and NNLO perturbation series and beta-function coefficients with, in the MS scheme, 

77139 - 15099iV/- + 325A^? 

= i728ir 

Since the NNLO r2 is unknown, so is p2- The coefficient Kq is directly related to the coefficient 
A of the 1/Q power corrections in eq. (^). 

Since R{Q) and dR/d\ogQ are both observables one could in principle directly fit eq. (|5|) 
to the data and thus constrain the unknown coefficients p2 and Kq. At asymptotic energies all 
observables run universally according to, 

'^^ -bR\l + cR), (9) 



dlogQ 



and one could see how close the data are to this evolution equation. Given the error bars of the 
data and the separation in Q of the different experiments it is preferable, however, to integrate 
up eq. (|5|) using asymptotic freedom {R{Q) — > as Q — > oo) as a boundary condition. In this 
way one obtains, 

1 + clogf^^^ =blog(-^] - r dx ( — + — -i , (10) 

where A^j is a constant of integration. By comparing with the Q — > oo behaviour of eq. @ one 
can deduce that A^j is related to A-jg^, 

where r = r^^{^ = Q). 

If we assume that the right-hand side of eq. (|5|) is adequately parameterized by, 

- bp{R) = -bR^{l + cR + p2R^) + KoR-^/^e-^^^^, (12) 

we can then insert this form into eq. ( |lO| ) and by (numerically) solving the transcendental 
equation, perform fits of p2, Kq (or equivalently A) and A;^ to the data R{Q). 

Fig. ^(a) shows the fit to the data (dashed line) obtained by setting p2 = X = 0. This 
corresponds to the universal running of the observable given in eq. (^). The fitted value is 
A^^ = 266 MeV with a x^/d.o.f = 81.7/32. We clearly see that the data is falling much too 
quickly with increasing Q for the asymptotic behaviour to have set in at these scales. The data 
favours a more steeply falling evolution which could be caused by either higher order corrections 
with a positive p2, or power corrections with non-zero Kq. We therefore perform a 3-parameter 
fit allowing p2, Kq and Aj^ to vary independently which is shown as a solid line in Fig. |2|(a). 
The minimum \^ fit is acceptable (x^/d.o.f = 40.4/30) and estimating an error by allowing 
within 1 of the minimum gives, 

= 2451^2 MeV with p2 = -16 T H and A = 0.27l|]-J^ GeV. 

These values of p2 and A are reasonably small, thereby lending support to our critical assumption 
that the evolution equation could be parameterized in this way. Converting the extracted value 
of A-jjg into as{Mz) we find, 

as{Mz) = 0.1194 ± 0.0014. (13) 

We see that our central value is remarkably close to that obtained by ElB. The main 
difference is in how the errors are determined. In our approach, the errors are estimated by 
allowing the uncalculated higher orders to be fitted by the data and the data prefers these to be 
small. In particular, the renormalization group scale-dependent logarithms are automatically 
resummed to all orders on integrating eq. (|5|) and do not add a spurious extra large uncertainty 
in the extraction of A-p^ (or equivalently as{Mz))- As higher order corrections become known, 
the new RS-invariant terms, p2^ ps etc., can be incorporated and the fit refined. 

With such an accurate value of as{Mz), we should expect that applying this approach to 
other variables should yield consistent results. Unfortunately, the method described here relies 
on having reliable data over a wide range of Q values. This limits its use to a very few variables, 
like thrust or the heavy jet mass. If we repeat the same analysis for the average of the heavy 
jet mass, <m'jj/s>, we find that a one parameter fit with p2 = X = gives a very poor fit, 
X^/d.o.f = 213/29 and A-p-^ = 368 MeV. As seen in Fig. |(b) the data evolves much faster 
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Figure 2: The average value of (a) 1-Thrust and (b) heavy jet mass obtained experimentally compared with the 
expectation of eq. (Rw- The dashed lines show the fit to the data with p2 — \ = while in (a) the result of the 
three parameter fit (to ^jfg, p2 and A) is shown as a solid line. In (b), the solid line shows the result of a two 
parameter fit (to p2 and A) using the value of A-j^g obtained from the three parameter fit to <f — T>. 



than the QCD prediction. However, ahowing both p2 and A to vary while using the value of 
~ MeV obtained from the <1— T> analysis gives a much more satisfactory description.^ 
Here, x^/d.o.f = 40.3/27 while p2 = 13 and A = 0.11 GeV are sufficiently small to support our 
choice of parameterization. 



4 Outlook 

The last decade has seen enormous progress both experimentally and theoretically in understand- 
ing high energy hadronic events. In many cases, we now have a quantitative understanding of 
the rates for strong interaction processes and detailed predictions for the structure of the events. 
One of the main theoretical uncertainties remains the renormalization scale uncertainty and the 
value of the strong coupling. New techniques for next-to-next-to-leading predictions are in sight 
and will help to reduce this error. However, a more immediate improvement may be obtained 
by resumming the ultraviolet logarithms explicitly. 
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